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We present the first complete optimization of quantum tomography, for states, POVMs, and 
various classes of transformations, for arbitrary prior ensemble and arbitrary representation, giving 
corresponding feasible experimental schemes. 
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A crucial issue in quantum information theory is the 
precise determination of states and processes. The pro- 
cedure by which this task can be accomplished is known 
as quantum tomography [H, @, 0|. 

The most general quantum measurement is described 
by a POVM, namely a collection of positive operators 
Pi G B(H) satisfying the normalization J2i^i = ^ 3- 
The probability distribution of the outcome i of the mea- 
surement is provided by the Born statistical formula 



Pi = Tr[pPi\. 



(1) 



Tomographing an unknown state p of a quantum system 
means performing a suitable POVM {Pi} such that every 
expectation value can be evaluated from the probability 
distribution pi = Ti[pPi]. In particular the expectation 
value of an operator A can be obtained when it is possible 
to expand A over the POVM as follows 



A = j2fdm, 



(2) 



fi[A] denoting suitable expansion coefficients. The ex- 
pectation of A is then obtained as (A) = J2i fi[A](Pi) ■ 
When expansion ^ holds for all operators BiTL) — i. e. 
B(Tt) = spanjPi} — the POVM is called informationally 
complete 3, Sl- 
it is convenient to associate every operator A G B(H) 
to a bipartite vector in Ti ® H in the following way 

d d 

A= A mn \m)(n\^\A)) = ^ A mn \m)\n). (3) 



Information-completeness of the POVM along with con- 
vergence of the series ^ rewrite as follows 
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\Af 2 <J2\((PAA))\ 2 <b\\A\\h A e B(H), (4) 



with < a < b < oo. Sets of vectors |P,)) satisfying 
condition Q are known as frames 0]. This condition 
is equivalent to invertibility of the frame operator F = 
Si The expansion in Eq. ^ can be written 

as follows 



in terms of a dual frame {Di}, namely a set of operators 
satisfying the identity ^2i\Pi))((Di\ = I. For linearly 
dependent frame {Pi} the dual {Di} is not unique. 

The request for the POVM {Pi} to be informationally 
complete can be relaxed if we have some prior informa- 
tion about the state p. If we know that the state belongs 
to a given subspace V C B(H) the expectation value is 



(A) = ((p\A)) = ((p\Q v \A)) 



(6) 



Qv orthogonal projector on V, whence the set {Pi} is 
required to span only V. 

For the estimation of the expectation (A) of an observ- 
able A, optimality means minimization of the cost func- 
tion given by the variance 6(A) of the random variable 
((Di\A)) with probability distribution Tr[/jPj], namely 



6(A) :=^|((A:|A))| 2 Tr[pP]-|Tr[M]p 



(7) 



In a Bayesian scheme the state p is randomly drawn from 
an ensemble S = {pk,Pk} of states pk with prior proba- 
bility pk, with the variance averaged over S, leading to 

Ss(A) := ]T |«A|^})| 2 Tr[p s Pi\ -£ Pjfc | T r [p k A}\ 2 (8) 



where ps = J2k PkPk- Moreover, a priori we can be inter- 
ested in some observables more than other ones, and this 
can be specified in terms of a weighted set of observables 
Q = {An, q n }, with weight q n > for the observable A n . 
Averaging over Q we have 



]T((A|G|A)) Tr[p s Pi\ - Y,Phq n \ Tr[p k A n ]\ 2 



(9) 

where G = J2 n 1n\A n )}((A n \. The weighted set Q yields a 
representation of the state, given in terms of the expec- 
tation values. The representation is faithful when {A n } 
is an operator frame, e. g. when it is made of the dyads 
corresponding to the matrix elements (j\p\i). 
Notice that only the first term of 8s,g depends on {Pi} 
and {Di}. If pi G V for all states pi € S, reminding Eq. 
((6|) the first term of Eq. |(9|) becomes 
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J2(( D i\QvGQv\Di}}Tr[psPi} 



(10) 



\A})=J2((Di\A)}\Pi}), 

i 



(5) 



We now generalize this approach to tomography of 
quantum operations, keeping generally different input 
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and output Hilbert spaces Hi n and H ou t, respectively. 
This has the advantage that the usual tomography of 
states comes as the special case of one-dimensional Hi n , 
whereas tomography of POVMs corresponds to one- 
dimensional H OU t- 

A quantum operation is a trace non increasing CP-map 
T : B(Tii n ) — ► B(H ou t)- In order to gather information 
about a quantum operation T, the most general proce- 
dure consists in: i) preparing a state p £ B(Hi n <S> Ha) 
where Ha is an ancillary system with the same dimen- 
sion of Hin\ it) measuring the state [T ® Ia){p) with 
a POVM {Pi}. The probability of obtaining a generic 
outcome i is given by 

Pi = Tr[(T®I A )(p)Pi\, (11) 

which, using the Choi-Jamiolkowski isomorphism [§], 

T(p) = Tr m [{I out ®p T )R T l R r = T®I in (\I))((I\) (12) 

becomes 

Tr[Tr m [(/4 <g> i?r)(A ® Iout)]Pi] = Tr^IT^], (13) 

where is the transposition w.r.t. the orthonormal basis 
in Eq. (J3J, and 

n[ p) = {Tr A [(p® I ut){hn ® P°° ut )]} T - (14) 

It is convenient to use here the notion of tester along with 
the theoretical framework introduced in [{J. A tester is 
the natural generalization of the concept of POVM from 
states to transformations, and is represented by a set of 
positive operators {IT} with 

Y^h = 1 ® °, Tr H = 1 (15) 

i 

The probability distribution in Eq. Jl3j is precisely 
represented by a Born- rule with the tester {IT} in place 
of {Pi}, and the operator Rr in place of p. Such gener- 
alized Born rule can be rewritten in terms of the usual 
one as follows 0| 

Pi = Tr[i? T IT] = Tr[T®l(u)Pi], (16) 

with 

v = |V5»«vH Pi = (Im-^Mma- 1 / 2 ). (17) 

This method allows a straightforward generalization of 
the tomographic method from states to transformation. 
Now tomographing a quantum operation means using a 
suitable tester IT such that the expectation value of any 
other possible measurement can be inferred by the proba- 
bility distribution pi = Tr[i?7-IT] . In order to achieve this 
task we have to require that {IT} is an operator frame 
for B(H ut <£> Hin)- This means that we can expand any 
operator on H ou t ® Hi n as follows 

A = ^2((A i \A))Il i AeB(H out ®H m ). (18) 



where {A^} is a possible dual of the frame {IT}, that is 
the condition £V |ir))((Ai| = I out ® I in holds. 

Optimizing the tomography of quantum operations 
means minimizing the statistical error in the determi- 
nation of the expectation of a generic operator A as in 
Eq. ((18)1 . This is provided by the variance 

8(A) = £ |((A^))| 2 Tr[i? r IT] - | Tt[R t A}\ 2 (19) 

i 

We assume an ensemble £ = {Rk,Pk} of possible trans- 
formations and a weighted set Q — {A n ,q n } of possible 
observables. Averaging the statistical error over these 
ensembles we obtain 

h,A ■= Y^{{Ai\G\Ai)) TrlRelli] - 5> fc g„| Tr[R k A n }\ 2 . 

i k,n 

(20) 

Optimizing this figure of merit means: i) optimizing the 
choice of the dual frame {A;}; ii) optimizing the choice of 
the frame {IT}. The optimization of the set {IT} reflects 
in both choosing the best input state for the quantum 
operation and the best final measurement. 

In the following, for the sake of clarity we will consider 
dim(Hi n ) = dim(Hout) =■ d, and focus on the "symmet- 
ric" case G = I; this happens for example when the set 
{A n } is an orthonormal basis, whose elements are equally 
weighted. Moreover, we assume that the averaged chan- 
nel of the ensemble £ is the maximally depolarizing chan- 
nel, whose Choi operator is Rg = g? -1 /® /. 

With these assumptions the relevant term of figure of 
merit becomes 

V = TrpT]. (21) 

i 

Since Rs is invariant under the action of SU(d) x SU(d) 
we now show that it is possible to impose the same co- 
variance also on the tester without increasing the value 
of r). Let us define 

m, g ,h ■■= (U g <8> V h )Ui{Ul 8 Vl), (22) 
&i, g ,h ■= {U g ® V h )Ai{Ul ® Vl). (23) 

It is easy to check that &-i, g ,h is a dual of Hi, g ,h by eval- 
uating the group average after the sum on i. Then we 
observe that the normalization of Tli.g.h gives 

Y^Jdgdh^g^^d- 1 !®! (24) 

i 

corresponding to a = d~ x I in Eq. (fTT)) . namely one can 
choose v = rf _1 |7))((/|. In the last identity dg and dh are 
invariant measures normalized to unit. 

It is easy to verify that the figure of merit for the 
covariant tester is the same as for the non covariant 
one, whence, w.l.o.g. we optimize the covariant tester. 



The condition that the covariant tester is information- 
ally complete w.r.t. the subspace of transformations to 
be tomographed will be verified after the optimization. 

We note that a generic covariant tester is obtained by 
Eq. 11221) . with operators IT becoming "seeds" of the co- 
variant POVM, and now being required to satisfy only 
the normalization condition 



E^w 



(25) 



(analogous of covariant POVM normalization in (3. fl(il|). 
The problem of optimization of the dual frame has been 
solved in With the optimal dual, the figure of merit 
simplifies as 



where 



^ = E Jdgdh 



Tr[n i; 



g,h\ 



(26) 



= ldgdhW g , h XWl h 



(27) 



with W q h 



U q ® U* <g> V h 



V£ and X 



d\Hi)) ((IT |/ Tr[IIj]. Using Schur's lemma we have [ij 



X : 

Pi 
P3 



Pi + AP 2 + BP 3 
: ^13 ® (I24 



CP 4 , (28) 
P 2 = {ha - flis) ® ^24, 

n 2i ) , p 4 - (/13 - fiia) ® (^24 - n M ), 



having posed O 
1 



|J))((/|/d and 



i4 = 
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d 2 - 1 
1 
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y, Tr[(Tr 2 jIT]) 2 ] 



E 



Tr[ry 

Tr[(Tn[n z ] 
Trpi] 

dTr[H|] 



- 1 



(29) 



(d 2 
One has 



1) 



E 



Tr[n ( 



(d 2 - 1)(A + B) - 1 



Tr^ 1 ] = 1 + (d 2 - 1) 



.4 



1 



(d 2 1) 
C 



(30) 



We note that if the ensemble of transformations is con- 
tained in a subspace V C B{H ou t ® Win) the figure of 
merit becomes r] = Ti[X$Q\>], where X* is the Moore- 
Penrose pseudoinverse. We now carry on the minimiza- 
tion for three relevant subspaces: 

Q = B(Hout <8> Wi„), C = {i? £ Q, Tr out [i?] = I m } 
U = {ReQ, Tr out [i?] = 7 in , Tr m [i?] = I out } (31) 

corresponding respectively to quantum operations, gen- 
eral channels and unital channels. The subspaces C and 



A, 
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Figure 1: Physical implementation of optimal quantum trans- 
formation tomography. The two measurements are Bell's 
measurements preceded by a random unitary. The state |^/)} 
depends on the prior ensemble. 



U are invariant under the action of the group {Wg,fe} and 
thus the respective projectors decompose as 



Qc = Pi+P2 + Pa, Qu = Pi+Pa 



(32) 



Without loss of generality we can assume the operators 
{IT} to be rank one. In fact, suppose that IT has rank 
higher than 1. Then it is possible to decompose it as 
II = with Ilj.j rank 1. The statistics of IT can 

be completely achieved by ITj through a suitable post- 
processing. For the purpose of optimization it is then 
not restrictive to consider rank one IT, namely IT = 
cti\ i S?i))(( i S?i\, with YiO-i = d. Notice that all multiple 
seeds of this form lead to testers satisfying Eq. ([25|) . 

In the three cases under examination, the figure of 
merit is then 



rt Q = Tr^ST 1 ] = 1 + (d 2 - 1 
Vc = Tr[XiQ c ] = l + (d 2 -l) 



2 {d 2 - 1) 



A 1-2A 

1 , (d 2 -l) 2 
A 1-2A 



f(d 2 -l) 2 

m = Tt[X$Qu] = l + {d 2 -l){ \_ 2 ^ 



where < A = {d 2 - l)" 1 ^ ^ Tr[(*itft) 2 ] 



(33) 



1) < 



< ^- The minimum can simply be determined by 
derivation with respect to A, obtaining A = l/(d 2 + 1) 
for quantum operations, A = l/(\/2(d 2 — 1) + 2) for gen- 
eral channels and A = for unital channels. The corre- 
sponding minimum for the figure of merit is 

r] Q > d 6 + d 4 - d 2 

nc >d 6 + (2V2 - 3)rf 4 + (5 - 4V2)rf 2 + 2(V2 - 1) 
m > (d 2 - l) 3 + 1. (34) 

The same result for quantum operations and for unital 
channels has been obtained in [13j| in a different frame- 
work. 
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These bounds are simply achieved by a single seed 
n = d| |, with 



Tr[(** T ) 



2d 



V2(d 2 



1) 



d 2 + r d(V2(d 2 -l) + 2)' 



1 (35) 



respectively for quantum operations, general channels 
and unital channels, namely with 



(36) 



where /3 = \J{d+ l)/(d 2 + 1) for quantum operations, 
P = [(d-l)(2 + V2(d 2 -1))]" 1/2 for general channels and 
(3 = for unital channels, and is any pure state. The 
informationally completeness is thus verified a posteriori 
(see 0). 

The same procedure can be carried on when the op- 
erator G has the more general form G = g\P\ + 52-P2 + 
ff3p3 + g^Pi, where Pi are the projectors defined in ((2Sj) . 
In this case Eq. ([30)1 becomes 



Tr[ir J G] = gi 



(d 2 - ljgj 
C 



(37) 

which can be minimized along the same lines previously 
followed. G has this form when optimizing measuring 
procedures of this kind: i) preparing an input state ran- 
domly drawn from the set {U g pU^}; ii) measuring an 
observable chosen from the set {UhAU^}. 

We now show how the optimal measurement can be 
experimentally implemented. Referring to Fig. Q], the 
bipartite system carrying the Choi operator of the trans- 
formation is indicated with the labels Si and S2. We 
prepare a pair of ancillary systems Ai and A2 in the 
joint state |^)}((^|, then we apply two random unitary 
transformations Ui and U% to Si and 52, finally we per- 
form a Bell measurement on the pair A\Si and another 
Bell measurement on the pair A2S2. This experimental 
scheme realizes the continuous measurement by random- 
izing among a continuous set of discrete POVM; this is a 
particular application of a general result proved in 
The scheme proposed is feasible using e. g . the Bell mea- 
surements experimentally realized in [15|. We note that 
choosing maximally entangled (as proposed for ex- 
ample in [16(]) is generally not optimal, except for the 
unital case. 

With the same derivation starting from Eq. (|2lf . but 
keeping dim(Hm) 7^ dim(H ou t) , one obtains the optimal 
tomography for general quantum operations. The special 
case of dim(H l „) = 1 (one has P 3 = P 4 = in Eq. ([28]) ) 
corresponds to optimal tomography of states, whereas 
case dim(H ut) = 1 (P-2 = P4 = 0) gives the optimal 



tomography of POVMs. The corresponding experimen- 
tal schemes are obtained by removing the upper/lower 
branch for POVMs/states, respectively. In the remaining 
branch the bipartite detector becomes a mono-partite, 
performing a von Neumann measurement for the qudit, 
preceded by a random unitary in SU(d). Moreover, for 
the case of POVM, the state |\&)) is missing, whereas, for 
state-tomography, both bipartite states are missing. The 
optimal 77 in Eq. (fT0[) is given by 77 = d 3 + d 2 — dAn both 
cases (for state-tomography compare with Ref. [Uj). 



In conclusion, we presented a general method for op- 
timizing quantum tomography, based on the new notion 
of tester. The method is very versatile, allowing to con- 
sider arbitrary prior ensemble and representation. We 
provided the optimal experimental schemes for tomogra- 
phy of states and various kinds of process tomography, 
giving the corresponding performance, all schemes being 
feasible with the current technology. 
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